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Abstract 

The dark photon, an new hypothetical light spin 1 field, constitutes a well-motivated dark matter candidate. It manifests as an 
oscillating electric field with a fixed direction, which can be observed in magnetometric records. In this letter, we use magnetometer 
data from the Voyager probes to look for the dark photon in the 10 24 eV to ICC 19 eV mass range, corresponding to frequencies 
between 10 9 Hz and 10 4 Hz. We also discuss the sensitivity of possible future SQUID magnetometry experiments. 
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1. Introduction 

We now have a compelling set of evidence for the existence 
of some form of cold, non-baryonic dark matter exerting grav¬ 
itational influence from astronomical to cosmological scales. 
Elucidating the microscopic nature of the dark matter is one of 
the most pressing challenges in physics. This important open 
question connects particle physics and cosmology. It is likely 
that future progress in one field will proceed from the advance 
in the other. Among the countless theories to explain the dark 
matter, two broad classes can be distinguished: WIMPs and 
WISPs. 

WIMPs are weakly interacting massive particles, with a typ¬ 
ical mass in the GeV to TeV range. They are still actively 
searched for by direct detection with low background detectors 
burned underground. Alternatively, the particle could be pro¬ 
duced at the LHC in high energy proton-proton collisions, it is 
looked for in events with large missing energy. 

WISPs are weakly interacting slim particles whose mass can 
be much less than 1 eV (see CO El for recent reviews). This 
type of dark matter is necessarily bosonic, since the Pauli ex¬ 
clusion principle forbids fermions to condense above a cer¬ 
tain density, the so-called Tremaine-Gunn bound 0. Possi¬ 
bly, these bosonic fields are produced in the early Universe via 
the vacuum misalignment mechanism 00. They behave very 
much like classical fields, oscillating at a frequency / = Mc 2 /H, 
where M is the mass of the bosonic field. Although the particles 
are very light as compared to the temperature of the early Uni¬ 
verse, the field behaves as cold dark matter. More precisely, the 
equation of state of the energy density stored in the oscillations 
is the same as that of pressureless dust. 

The scenario of dark matter as a condensate of very light 
bosons has mostly been discussed in the framework of the QCD 
axion, a spin 0 field invented to solve the strong CP problem 0 
|5]. Theory predicts that dark matter axions could be converted 
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into microwave photons in the presence of an external magnetic 
field. Haloscopes, i.e. devices realizing this conversion such 
as ADMX 0, are probing a large portion of plausible axion 
masses and couplings. 

Recently, another form of WISPy dark matter has been pro¬ 
posed by Nelson and Scholtz (7J. The scenario consists in 
adding a new vector field associated with a hidden 1/(1) gauge 
symmetry called the dark photon. The energy stored in the os¬ 
cillations of this field would constitute the dark matter. In this 
letter we explore the possibility of detecting these oscillations 
using precise magnetometers. The proposed technique could 
be applied in a broad frequency domain ranging from 10 9 Hz 
to 1 kHz, corresponding to dark photon masses in the range 
10 24 - 10 12 eV. It is thus complementary to the newly pro¬ 
posed search with a dish antenna in the radio domain ®. The 
paper is organised as follows. In section [2] we review the dark 
photon dark matter scenario and derive the expected magnetic 
signal. In section [3] we set a bound on the kinetic mixing x 
using the magnetometric data recorded by the two Voyager in¬ 
terplanetary missions. Then, the sensitivity of a dedicated mea¬ 
surement with SQUID magnetometers is discussed in section 

m 

2. Magnetic signature of dark photon dark matter 

We consider an extension of standard electrodynamics with 
an additional spin 1 field X fl 0, defined by the Lagrangian den¬ 
sity 

1 1 M 2 

Z = ~-F IJV F^' - -X MV X» V + — X^ (1) 

where F )n , = d M A v - d v A t , represents the field strength of the 
ordinary photon field ,4 /( , J l ‘ is the usual current of charged 
fermions, and X flv = d t ,Xv - is the field strength of the new 
boson of mass M. The name of the new boson is not unique 
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in the litterature, it was first referred to as a paraphoton, then 
more recently as hidden photon or dark photon. Since we are 
regarding the new field as a dark matter candidate we choose to 
use the latter name. 

The dark photon couples to the rest of the world through the 
kinetic mixing term with the dimensionless coupling^. String- 
inspired extensions of the Standard Model of particle physics 
typically predict a value for x in the range 10 12 - 10 3 IfTOl . 
By making the transformation ,4 /( = A M + xX^ it is apparent 
that a particle of charge e for the usual photon has a suppressed 
effective charge of x e f° r the dark photon. 

The classical equations of motion for the dark photon field 
admit oscillating solutions of the form 

X i = Re(x i e- iM ‘) , X Q = 0, (2) 

where x, are complex amplitudes. The corresponding complex 
amplitude for the dark electric field -8qX is £dark = iMx. The 
energy density stored in the oscillation is p = ^-\x\ 2 . In terms 
of the electric field, in SI units we have 

P = ll^darkl 2 - (3) 

Following Nelson and Scholtz S3 we attribute the energy 
density of dark matter to the oscillations of the dark photon. 
We use the canonical value of p = 0.3 GeV/cm 3 ED for the 
local dark matter density and find £dark = 3300 V/m. 

This dark electric field oscillates in the rest frame of the dark 
matter. Our motion relative to the halo induces a motional dark 
magnetic field with complex amplitude given by 

v -> 

Bdark = T X B'darki (4) 

C z 

where v is the relative velocity of the Earth in the halo, that 
we assume to point to the Cygnus constellation. We use the 
canonical value for solar velocity v = 254 km/s. 

The dark magnetic field couples to the usual electric charges 
with the coupling constant scaled by e — * x e - In other terms, 
we should search for a magnetic field B = xB&a, oscillating at 
the frequency / = Mc 2 /h, in a plane transverse to v, with an 
amplitude of 

B - x sin(0)v£ dar k/t' 2 = 9.3 x sin ^) nT - ( 5 ) 

where 9 is the angle between our velocity v and the direction of 
the dark electric field. 

3. Constraint from the Voyager magnetic probe 

In order to search for the signal described above in the low- 
frequency regime (about 10 9 Hz or a mass of around 10 22 
eV), we looked for long-running magnetic experiments in low 
noise environments. One such experiment is the Voyager Pro¬ 
gram. The two Voyager probes were constructed by the Jet 
Propulsion Laboratory and launched in 1977. Among many 
other instruments, they carried triaxial fluxgate magnetometers 
El. The last planet encountered by Voyager 1 was Saturn in 





Figure 1: Magnetic field values as a function of time. The colored parts are 
the the data actually used, running from 1990 to 2011. The blue line is the 
data from Voyager 1, the red one from Voyager 2. The upper figure show the 
magnetic field along the x (apical) axis, the middle figure along the y axis and 
the bottom figure along the z axis. For readability, only one in ten points is 
actually shown. 


1980, whereas Voyager 2’s last flyby was above Neptune in 
1989. Between these events and their arrival at the termina¬ 
tion shock, respectively in 2004 for Voyager 1 and in 2007 for 
Voyager 2, the probes were deep in interplanetary space. 

The data set running from 1990 to 2011 presented a noise be¬ 
low 1 nT and was very well suited to our needs. The magnetic 
field seen by both probes during their entire flight is presented 
in figure [T] From now on, we will exclusively use magnetic 
data taken between 1990 and 2011. The data were taken from 
the Goddard Space Flight Center’s OmniWeb and FTP services. 
It consists of triaxial hourly samples of the external magnetic 
field. These values are given in RTN coordinates (Radial Tan¬ 
gential Normal); it is then necessary to transform them into a 
referential more adapted to our aims. There is clearly a pre¬ 
ferred direction which is along the speed of the probes in the 
galactic referential; this direction is essentially that of the speed 
of the Sun in the galaxy, also known as the solar apex. Along 
this axis, there can be no signal from the dark photon, and thus 
the corresponding magnetometric data component can be used 
to discriminate oscillations resulting from sources other than 
the dark photon. This direction is called the apical direction or 
x axis, the plane orthogonal to it is called the transverse plane. 
We have then selected two other axes in order to get a direct tri¬ 
hedron; those are arbitrary in the sense that our results will not 
depend on this choice. We call them the y and z axes. For our 
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Figure 2: Noise spectral density of the transverse magnetic field (mean data 
of the two probes). The black line is the signal proper, the red line is a moving 
average calculated on 100 points. 

purposes, we have chosen a y axis in the galactic plane, pointing 
in the direction of positive galactic longitude and orthogonal to 
the solar apex, and a z axis completing the orthonormal coor¬ 
dinate system. In order to use the data of both probes, we will 
work with the mean of the magnetic values of Voyager 1 and 
Voyager 2 for each hour and each axis. The magnetic values in 
the tranverse plane are coded into a complex number of form 
s{n) = B ' yA ^ ly + 2 -' ; where B\ y is the component of the field 

along the y axis seen by Voyager 1 at sample time number n, 
and other values are labeled accordingly. We then compute the 
discrete Fourier transform of this signal, given by : 

N-l 

S(k) = Y J s(n)e~ 2i ^, ( 6 ) 

o 

where S is the fourier transform, s is the signal described above, 
and N is the number of points considered in the transform. Once 
we have this result, it is easy to get the noise spectral density of 

the signal, given by S x (k) = yj IW)* Af l~ ~, w j t j 1 /y t (jj e sampling 
period and T the total duration of the acquisition. Since the 
noise density depends only on the module of the Fourier trans¬ 
form, it is independant of the explicit choice of the y and z axes. 
The results are given in the figure [2] 

It is immediately apparent that there is a one day period 
(1.1 x 10 5 Hz) signal in the data, along with associated har¬ 
monics. This is by no means a sign of the dark photon, as the 
same periodicity appears in the apical field. This peculiar oscil¬ 
lation stems from the values missing in the record. The figure 
[^presents a simulated case with a 10 21 eV and^- = 1 dark 
photon, with distorsion induced by the points missing on the 
record. The expected mass peak at 2.4 x 10~ 7 Hz is present, but 
we see 24-hour components and associated harmonics arising 
as well. Those frequencies remain even if we get rid of any sig¬ 
nal and simply fill any non-missing data with a constant value. 
Thus, we can safely assume that those peaks are just an artifact. 
We can observe that other than these previously mentioned fre¬ 
quencies, no significant peak can be seen. 
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Figure 3: Noise spectral density of a simulated dark photon signal with param¬ 
eters M andx and distorsion induced by the missing data. 


Table 1: Constraints (at 95 % C. L.) on the amplitude of an oscillating magnetic 
field in apical and transverse directions in the frequency range accessible with 
Voyager magnetic data. 


Frequency 

Bj 

Ba 

lO” 8 Hz 

16 pT 

9 pT 

1(T 7 Hz 

11 pT 

7 pT 

1(T 6 Hz 

4 pT 

2.2 pT 

1(T 5 Hz 

1.5 pT 

0.7 pT 

1(T 4 Hz 

0.3 pT 

0.2 pT 


The next step is then to determine an exclusion limit of the 
values of ^fsin(0). Since the noise level varies with the fre¬ 
quency, no single criterion of exclusion can be used on the 
integrality of the bandwidth. Instead, we compute a moving 
average of the noise spectral density by blocks of 200 points 
(corresponding to a bandwidth of 3.0 x 10 7 Hz), along with 
the standard deviation of these blocks. Having verified that the 
noise is gaussian in nature, with 95 % confidence, we exclude 
any signal with a given frequency and an amplitude larger than 
the mean at this frequency plus two standard deviations. For 
low frequencies, this method is inadapted, as the noise levels 
vary greatly from one point to the next. For the 200 first points 
(from 1.5 x 10~ 9 Hz to 3.0 x 10~ 7 Hz), we exclude any signal 
greater than the noise spectral density of the data at a given fre¬ 
quency plus two times the standard deviation of these 200 first 
points. In order to do that, we simulate a dark photon oscillating 
signal with an arbitrary frequency and a value of the coupling 
parameter^ set to 1, and we compute the noise spectral density, 
which has a maximum value of 2.12 nT/VHz. This density 
does not depend on the mass of the particle to first order, and 
is linear in which allows us to transform the previous limit 
on the noise density into a limit on^sin(0). The final limit is 
represented on figure[4] compared to expected limits from other 
experiments. For completeness, the same analysis has been ex¬ 
tended to the apical field; results are given in table[T|in the form 
of constraints on the effective field amplitudes. 
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Figure 4: Presentation of exclusion limits on x as a function of the mass of the 
dark photon. Black solid line : this work. Red dotted line : expected constraints 
from a combination of radio observations ED. Blue solid line : Constraints 
from Cosmic Microwave Background observations 03 Black dotted line : 
SQUID magnetometry expected limit, explained in section[4] 
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4. Conclusion and prospects 

We have shown that precision magnetometry could be used to 
search for the dark photon, a newly proposed form of dark mat¬ 
ter. In this scenario, the dark matter is attributed to the coherent 
oscillation of the new held, which acts as a pseudo magnetic 
held due to our motion in the halo. The very long magnetic 
records from the Voyager 1 and 2 space probes were used to set 
competitive limits in the mass range 10 23 eV to 10 18 eV, cor¬ 
responding to frequencies between about 10 9 Hz and 10 4 Hz. 
Interestingly, this covers the scenario known as Fuzzy dark mat¬ 
ter 1131 . With 95 % confidence, we can exclude values of the 
coupling parameter^ as low as 3 x 10 s for the highest masses, 
which represents the best limit for this range of frequencies. 

In the future, a natural development would be to look at 
somewhat shorter periods. It is sensible to propose exploration 
of such periods (from the year to the millisecond) with a simi¬ 
lar approach, using SQUID (Superconducting Quantum Inter¬ 
ference Device) magnetometry. Using noise data published 
in Henry et al. ED, which gives a noise spectral density of 
4.0 fT/VHz above 1 Hz and of (4.9//) fT/Vfe, we suppose a 
data acquisition for a time of one month and that we can see a 
signal if its noise spectral density is two times as large as that 
of the ambient noise. We can then expect a sensitivity as low 
as x = 10 s • The predicted limit is entitled ’’SQUID magne¬ 
tometry” in figure [4] More generally, any precision magnetom¬ 
etry experiment, for example a clock comparison experiment 
D3, could be adapted to search for these oscillations. These 
prospects open new perspectives in the hunt for the dark matter. 
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